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1. Introduction and summary 

The AdS/CFT correspondence [1] has been the subject of intensive research. This duality 
provides a powerful tool to study a plethora of interesting topics in theoretical physics. 
One of the best studied examples is = 4 SYM gauge theory which is conjectured to be 
dual to superstring theory on AdSs x S^. However, even testing or proving the duality for 
these concrete models is a difficult problem. 

A breakthrough in the understanding of this duality was the discovery of integrable 
structures. Integrability was found in = 4 super Yang-Mills theory by the appearance 
of (integrable) spin chains describing the operator spectrum [2]. The classical string sigma 
model on AdSs x S^ was also shown to be integrable [3]. If the full quantum theories also 
exhibit integrable structures, then that would constrain them considerably. For example, 
in scattering processes for integrable theories the set of particle momenta is conserved 
and scattering processes always factorize into a sequence of two-body interactions. This 
implies that in such theories the scattering information is encoded in the two-body S- 
matrix. Unfortunately, a full proof of integrability of both theories is currently still lacking, 
but nevertheless there is a lot of evidence hinting that integrability is a feature of the full 
quantum theories. 

By assuming integrability one can make use of the S-matrix approach. This proved to 
be a powerful instrument to study the operator spectrum of AA = 4 SYM. [2, 4-6] . This lead 
to the conjecture of the "all-loop" Bethe equations describing the gauge theory asymptotic 
spectrum [7,8]. For the AdSs x S^ superstring, based on the knowledge of the classical 



- 1 - 



finite-gap solutions [9], a Bethe ansatz for the su(2) sector was proposed [10]. Finally, in 
both cases, exact two body S-matrices were found [8, 11]. These enabled the use of the 
Bethe ansatz [8, 12, 13], confirming the conjectured Bethe equations for physical states. 

More precisely, the two body S-matrix (scattering fundamental multiplets) that ap- 
pears in this approach, is almost completely fixed by symmetry. Both the asymptotic 
spectrum of A/" = 4 super Yang-Mills theory [5, 8] as well as the light-cone Hamiltonian 
[14-16] for the AdSs x supcrstring exhibit the same symmetry algebra; centrally extended 
su(2|2). The requirement that the S-matrix is invariant under this algebra determines it 
uniquely up to an overall phase factor [8] and the choice of representation basis [11]. In 
a suitable local scattering basis, the S-matrix respects most properties of massive two- 
dimensional integrable field theory, like unitarity, crossing symmetry and the Yang-Baxter 
equation [11]. Furthermore, the (dressing) phase appeared to be a striking feature of the 
string S-matrix and it has been studied intensely, see e.g. [17—20]. By combining its 
expansion in terms of local conserved charges with crossing symmetry [21], one can find 
interesting solutions [17, 22], which incorporates string and gauge theory data. 

Actually, apart from the multiplet of fundamental particles, the string sigma model 
also contains an infinite number of bound states [23]. These bound states fall into short 
(atypical) symmetric representations of the centrally extended 5u(2]2) algebra [24-26]. Of 
course, these states scatter via their own S-matrices. Recently a number of these bound 
states S-matrices have been found; they describe scattering processes involving fundamental 
multiplets and two-particle bound state multiplets [27] and the scattering of a fundamental 
multiplet with an arbitrary bound state multiplet [28]. However, invariance under centrally 
extended su(2]2) is not enough to fix all these S-matrices. One needs to impose the Yang- 
Baxter equation by hand in order to fix them up to a phase. The overall phase that still 
remains can be chosen to satisfy crossing symmetry [27]. 

Both the fundamental and the two particle bound state S-matrices were shown to 
exhibit a larger symmetry algebra ^ of Yangian type [29-32]. Moreover, as an alternative 
to the Yang-Baxter equation, the bound state S-matrices appear to be completely fixed 
(again up to an overall scale), by requiring invariance under Yangian symmetry [32]. This 
seems to indicate that the Yangian symmetry is restrictive enough to fix all bound state 
S-matrices and perhaps should be seen as the fundamental scattering symmetry. 

The formulated asymptotic Bethe ansatze, following from the S-matrix approach, only 
describe the spectra in the infinite volume limit. However, for a full check of this par- 
ticular case of the AdS/CFT correspondence, the complete spectra have to be computed 
and compared. Away from the asymptotic region, matters become more involved because 
wrapping interactions appear. One way to include these is Liischer's perturbative approach 
[33-36]. In Liischer's approach one deals with corrections coming from virtual particles that 
propagate around the compact direction. These virtual particles can be both fundamen- 
tal particles and bound states. This method has proven quite successful since it recently 
allowed the computation of the full four-loop Konishi operator, including wrapping inter- 

^See [11] for an earlier discussion of higher symmetries of the fundamental S-matrix. 
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actions [28, 37]. The result coincided with the gauge theory computation [38-40], providing 
an extremely non-trivial check of the correspondence. 

The approach by Liischer is closely related to the thermodynamic Bethe ansatz (TBA) 
[41]. In the TBA one deals with finite size effects by defining a mirror model [23]. Finite 
size effects in the original theory correspond to finite temperature effects in the infinite 
volume for the mirror theory. Here again, one needs to include all (physical) bound states. 
One of the advantages of this method is that one can still use the asymptotic Bethe ansatz. 
So far, this has only been carried out for the fundamental multiplets [8, 12, 13]. 

In other words, knowledge of the bound states, their S-matriccs and the corresponding 
Bethe ansatze is crucial for a complete understanding of finite size effects. Because of the 
relation between Yangian symmetry and the bound state S-matrices, one might suspect 
that the Bethe equations are also closely related to Yangian symmetry. This indeed appears 
to be the case as we will explain in this note. The Bethe equations that are derived can be 
used to study bound states of the mirror model by analytic continuation [23]. 

The aim of this paper is to provide a rigorous derivation of the Bethe Ansatz equations 
for the bound states, by diagonalizing the multi-particle bound state S-matrices using 
Yangian symmetry as a tool. The explicit bound state number dependence in the Bethe 
equations appears through the parameters and the dressing phase only. The Bethe 
equations emerging in this way, were found, a posteriori, to follow from a fusion procedure. 
This justifies the fusion procedure at the level of the Bethe ansatz equations. 

The paper is organized as follows. In the first section we recall the basic formulas 
of the Yangian of the centrally extended 5u(2|2) algebra in the superspace formalism. 
Subsequently, we will discuss the Bethe ansatz for fundamental representations. Finally, 
the Bethe ansatz will be reformulated in terms of coproducts of the Yangian and the bound 
state Bethe equations will be derived. 

2. Symmetry algebra, coproducts and S-matrices 

In this paper the full (Yangian) symmetry algebra of the S-matrix is important. In this 
section we will give a brief overview of the Yangian symmetry of the bound state S-matrices. 
For more details see e.g [32] and references therein. 

2.1 The algebra in superspace 

We will first discuss centrally extended 5u(2|2). This algebra is the symmetry algebra of the 
AdSs X S^ superstring and it is also the symmetry algebra of the spin chain connected to 
A/" = 4 SYM. The algebra has bosonic generators M, L, supersymmetry generators Q, G and 
central charges H, C, C^. The non-trivial commutation relations between the generators are 
given by 





(2.1) 
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The eigenvalues of the central charges are denoted by C, . The charge H is Hermitian 
and the charges C, and the generators Q, G are conjugate to each other. 

For computational purposes, it proves worthwhile to consider representations of the 
algebra in the superspace formalism. Consider the vector space of analytic functions of two 
bosonic variables wi^2 and two fermionic variables ^3^4. Since we are dealing with analytic 
functions we can expand any such function ^{w, 9): 

00 

r'-'''-'''^Wa,...Wa,_,e^e0. (2.2) 

The representation that describes ^-particle bound states of the light-cone string theory on 
AdSs X of centrally extended su(2|2) is 4^ dimensional. It is realized on a graded vector 
space with basis |e(jj...a^), |cai...a^_ia ), |eai...a^_2a/3), where are bosonic indices and a, {3 are 
fermionic indices and each of the basis vectors is totally symmetric in the bosonic indices 
and anti-symmetric in the fermionic indices [24, 25, 27]. In terms of the above analytic func- 
tions, the basis vectors of the totally symmetric representation can evidently be identified 
|eai...a«) ^ . . . Wai,\ea^,„a^_^a) Wa^ . . . Wa^_-^6a and |eai...a<._ia/3) ^ Wai • • • UJae^i^aOp 

respectively. In other words, we find the atypical totally symmetric representation de- 
scribing ^-particle bound states when wc restrict to terms For later convenience, we 

introduce the notation W^'^]r^^.^i = {w^^^ y- {w^^ y {9^^'^^ )^ (9^^^ Y , where (m) denotes differ- 
ent copies of the representation. Clearly, for an ^-particle bound state representation, this 
means that i+j + k + l = i and k,l = 0, 1. 

In this representation the algebra generators can be written in differential operator 
form in the following way 



f b d _ Ixb,,, d Ta — a d _ IxPa d 

and the central charges are 

n = {ad + hc) L,^^+6„4-^ 



(2.3) 



(2.4) 



To form a representation, the parameters a,b,c,d satisfy the condition ad — be = 1. The 
central charges become i dependent: 

H = £{ad + bc), C = £ab, = £cd. (2.5) 

The parameters a, b, c, d can be expressed in terms of the particle momentum p and the 
coupling g: 

9_ V ^_ /~g" x+ 

2£ CFF " ~ Y 2£ jr? 
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where the parameters satisfy 

x+ + ^-x-- — = —, — = e'P (2.7) 
x+ X g X 

and the parameters rj are given by 

r) = e'^r]{p), ■q{p) = e^^^/ix- - ix+, ( = e^'^. (2.8) 

The fundamental representation, which is used in the derivation of the S-matrix scattering 
fundamental multiplets [8, 11] is obtained by taking £ = 1. 

2.2 Yangians, coproducts and the S-matrix 

The double Yangian DY{q) of a (simple) Lie algebra g is a deformation of the universal 
enveloping algebra U{g[u,u~^]) of the loop algebra gfujU"^]. The Yangian is generated by 
level n generators Jj^, n G Z that satisfy the commutation relations 

KJn]=FS^SZ+n + 0{h), (2.9) 

where F^^ are the structure constants of g. The level-0 generators JJ^ span the Lie-algebra. 
The su(2|2) Yangian can be supplied with a coproduct in the following way [31,42,43] 

. n— 1 

n 



m=0 

A{U)=U0U, (2.10) 



where U comprises a so-called braiding factor. 

An important representation of a Yangian is the evaluation representation. This rep- 
resentation consists of states \u), with action S^W) — u"'Sq\u). In this representation the 
coproduct structure is fixed in terms of the coproducts of Jo, Ji- For the remainder of this 
paper we will work in this representation and identify Ji = J = ^ttj for the su(2|2) Yangian. 
One also finds that u is actually dependent on the parameters x"^ via Uj = xf + ^ — 

The S-matrix is a map between the following representations: 

S : YeM, (^'n ® -^iM, 1) nAPi, 1) ® '^h(P2, e'n, (2-11) 

where '5^i(pi, e^'^) is the ^^-bound state representation with parameters ai,bi,Ci,di with 
the explicit choice of ( = e^'^. This specific choice removes the braid factor U from the 
discussion [11]. 

The bound state S-matrices are now fixed, up to a overall phase, by requiring invariance 
under the coproducts of the (Yangian) symmetry generators 

S A(J^) = A°f (J^) S 

S A(J^) = A°P0^) S, (2.12) 
where A°p = PA with P the graded permutation. 
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For completeness and future reference, we will give the explicit formulas of the different 
coproducts and of the parameters Oj, fej, Cj, dj. First, the su(2|2) operators: 

A(J^)=J^o+4o- (2-13) 
The coproducts of the Yangian generators are 

- \^'0^i-,^^,c (2-14) 

and for the central charges 



A(M) 
A(C) 
A(Ct) 



4 



4^ 



'1^2 
1„ 



Ml + H2 + CiC^ 

Ci + C2 + ^HiC2- 



'1^2 



(2.15) 



The product is ordered, e.g. Q1Q2 means first apply Q2, then Qi (as differential operators). 
Finally, the coefficients used in A are given by: 



ai = 




Cl = 








m = 


e* 4 e 2 



0.2 
C2 
'72 



9 VI 
2^1 x+ 



^ ixi — ixf 



2f^2 

9 m 

/■ - 

e 4 4 / tXr, 



h 

di 
d2 



V 2h y^- 

I 9 ^2 ( "^2 1 I 

V V4 ~ ) 



(2.16) 
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The coefficients in A"*' are given by: 




The non-trivial braiding factors are all hidden in the parameters of the four representations 
involved. 



3. The su(2|2) coordinate Bethe Ansatz 

In this section we will briefly discuss the coordinate Bethe Ansatz for the su(2|2) string 
S-matrix as done in [8, 13]. 

3.1 Formalism 

Consider excitations with momentum pi, . . . , pj^i . We divide our space into regions 
V\Q, with V, Q permutations. In the sector where the particle coordinates are ordered as 
XQj^ < . . . < XQ^^ the ansatz for the wave function is given by 

=E / ci^{<l.^a,,e^^w}0«i(xi)...</,"-(x;,i) (3.1) 
■p 

The indices ai denote the type of the ith particle and </>"'(a:j) is a creation operator that 
creates a particle of type ai at position Xi. In general when the system contains bound 
states, the indices also run over bound states. The different sectors are related to each 
other by permutations and S-matrices. To be more precise 

= Si,,-A^'IS', (3.2) 

where P'\Q' are the permutations obtained from P\Q by interchanging the neighboring ith 
and _7th particles. To make this transformation property more explicit we will write 

<'.^a.^<'.^a.^ai---^a„ (3.3) 

where we do not sum over repeated indices. Periodicity of the total wave function is then 
formulated in the following Bethe equations 

Sfc^^lS := . . . S,^iS,i . . . = e^f'=^^^l2. (3.4) 
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The Bethe ansatz now consists of making a ansatz for the coefficients A^^^ in such a way 
that they solve ( |3.2| ). This ansatz is then plugged in the above formula and the explicit 
Bethe equations can be read off. 

Obviously, since any scattering process reduces to a product of two body scattering 
processes, we only need to consider two sites of the wave function. In what follows we will 
solve the for the coefficients A^^^- 

* 1 ■ ■ "'j^I 

3.2 Solving the coefRcients 

We will now solve the coefficients for the su(2|2) S-matrix for fundamental representations. 
In order make the derivation more transparent, we explicitly identify the coefficients A^^'^ 
with states |-) G '^i{pi) ® . . . '^(p/fi)- Fi'om now on we will omit the explicit mentioning 
of the sector "PIQ we are in. The different sectors are related via relabelling of momenta 
and interchanging of particle positions and hence the coefficients should be well-defined in 



the sense that they should respect this. From (3.2), we see that in order to be well-defined, 
we must have that 

m = 10)-, (3.5) 

where {(p)-,^ is the coefficient as constructed in the sector V'lQ'. 
Let us first define the vacuum: 

\0) =wi^K..wi^\ (3.6) 

One can also take other vacua but this leads to equivalent equations. If we, for the moment, 
consider the undressed S-matrix, we find that 

S|0) = |0). (3.7) 

In the light of ( |3.5D , we indeed find that this vacuum corresponds to a vacuum in all sectors 
V\Q and is well defined. 

The complete dressing phase So{pi,pj) can be easily incorporated by a rescaling 

|0) = 7 " • (3.8) 

'Ui^jSoiPuPj) 

This works because So{pj,pi) = So{pj,pi)^^ . To avoid making the formulas too cumber- 
some, we will restrict to the undressed S-matrix in the remainder of the paper. We will, of 
course, include the dressing phase in the full Bethe equations. 

The next thing to consider is the case where we have a fermion in this vacuum (the 
case in which the other boson is inserted is treated later on). We make an ansatz of the 
following form: 

\a) := My)wi'^ . . . . . . W^f ), ^k{y) = f{y,Pk) U S'''\y,pi). (3.9) 

i Kk 
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Figure 1: Example of the construction of the coefficients in the different sectors. The black dots 
correspond to the bosons of the vacuum, Wi, and the red dots correspond to the inserted fermions. 
The construction is well-defined if, by acting on the coefficient in the upper left quadrant by the 
S-matrix, one obtains the coefficient in the lower right quadrant. 



We will denote S{y,p) = S^^'^{y,p). We must check whether this construction is well 
defined in the sense that it respects (|3.5| ). How this works is schematically depicted in 
Figure 

By the factorization property of the S-matrix, it suffices to restrict to a two particle 
state and a two particle S-matrix. For this case (3.5) gives the following equations 



+ 2 ^ f{y,P2)S{y,pi) = f{y,P2) 



e 2 



viPi) xf — X2 ' xf — X 

^iE2.X^-X^ 'n{pi)x^-X2 



e--^ ^ + i. f{y,p2)S{y,p,) = f{y,p,)S{y,p2).{3.W) 

X]^ - X2 V\P2) Xj — X2 

These equations can be solved explicitly and the solution is given by: 

f{y,Pk) = v(.Pk)] 

S{y,Pk) = ih^^^-^, (3.11) 

where y enters as a integration constant. With a modest amount of foresight, we choose 
the overall normalization of / to be dependent on the bound state number ik- 

Actually this solution is the unique solution of ( p.5| ) for one fermion in the vacuum up 
to overall normalization (and hence the only well-defined coefficient for this case). This 
can easily be seen from ( p.lO| ). The general wave function is of the form 




a) 



AiTyW . . . + A^wl'^wi^') . . . Wf'^ + .... (3.12) 



Since we are only interested in solving ( |3.4D , the normalization of the state is irrelevant. 
Let us therefore pick a convenient normalization, Ai = f{pi), then by the first equation of 



( 3.10 ) one finds A2 in terms of Ai. This indeed fixes A2 as f{p2)S{pi) and by induction the 
rest of the coefficients also follow. Finally, let us stress that by fixing the normalization, 
one automatically finds from the S-matrix that the other coefficients of the wave function 
are of a f actor ized form. 
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The problem becomes more involved upon inserting two excitations. In this case, it 
appears that the ansatz for the coefficient can be written as 



|a/3) = |a/3)yiy2 + S"|a/3)j,iy2, (3.13) 
where we introduce a new S-matrix S^^ 

S"l«/3)?;iy2 = M[yi,y2)\aP)y^y^ + N{yi,y2)\Pa)y^y^. (3.14) 
Explicitly, it is given by 



\al3) = J2 ^fc(yi)^i(y2)H^i^'^ . . . VF« . . . w^'^ ...wi 



+ 



k<l 



^fc(yi)*z(y2)H^r^ ...w^^y.. w^'^ . . . lyf ) + (3.15) 



k<l 



+e-^ ^ ^k{yi)^k{y2)h{yi,y2,Pk)wi'^ . . . W^'^ . . . W^""^ 

k 

When we restrict to just two sites, the wave function splits into the sum of a wave function 
with just one fermion. By construction, this piece is exactly as described above. The other 
piece contains two excitations and is given by 

= {f{yi,Pi)f{y2,P2)S{y2,Pi)+Mf{y2,pi)f{yi,P2)S{yi,pi)]Wl^^WP 
+Nf{y2,Pi)f{yi,P2)S{yi,pi)wf^wl^^ 

+e"^/i(yi,y2,Pi)/(y2,Pi)/(yi,Pi)Tyi'Vf^ (3.16) 

+e-^/i(yi,2/2,P2)/(y2,P2)/(2/i,P2)5(y2,Pi)5(yi,Pi)Wi^'Vf 



Plugging this into ( |3.5| ) again allows one to find the explicit (unique) solutions of the 
unknown functions: 



M{yi,y2) 



2i/g 



y^ j/i 1/2 g 



N{yi,y2) = ^;t'' '\'\^ (3-17) 



yi + y-^_ y2 y^ g 



Kyi^y2,Pk) 



yiy2-xlx^ x^-xj^ yi - 2/2 



hr]{pkf ym x^ yi + ^ - y2 - ^ - f 



In this process, we introduced a new S-matrix and we also insist on naturalness of the wave 
function with respect to this S-matrix, S^^. One can repeat the above procedure to deal 
with this. For details we refer to [8, 13]. In this process we are again lead to introduce new 
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functions S^^'^^ , f^'^\ S^^^'^^ , S^^^'^^^ . The result of this consideration is 



S 



iiijr 



w,yk) 



-M -N 


= 1 






W -Vk- 


i 

a 


W -Vk + 


i 
9 


W -Vk- 


i_ 
9 


Wl — 102 - 


_ 2i 
9 



(3.18) 



The Bethe equations ( p. 41) are formulated in terms of the factors used in the ansatz in the 
following way 



III 

B=l 1=1 



{B,l)^{A,k) 



(3.19) 



where A, B denote the different levels and S^^ can be seen as the S-matrix describing 
scattering at different levels. Moreover, e*'^i''= = e~^^^'' and e*'^"''= = e^^™'^ = 1 are phases. 

By putting all of this together one obtains the well-known Bethe equations describing 
the asymptotic spectrum of the AdSs x superstring: 



n 



So{pk,Pi) 



^ 2 ^^(") _ (a) 

- Vl 



nn 



a=l 1=1 ^k 



(a) 




n 



Vk 



X, X 



(a) 

1=1 Vk ~ 




n 



Vk 



k 1=1 Vk ^ (a) 



w 



(a) 



(3.20) 



Kj\ „,,(") 1 I i K™ (a) (a) oi 

^k -Vk -^ + -9 wi'^^ - u;}") - f 

11 (a) _ (a) 1 11 



yi_ 

(q) (a) 1 
1=1 '^k Vk ~W ~ i l^k 



where a = 1,2 reflect the two copies of su(2|2) and So{pk,Pi) is the overall phase of the 
S-matrix. 

Notice that in this construction, the explicit form of the S-matrix is used. However, 
since not all bound state S-matrices are known, one must approach this problem in a 
different way if one wishes to find their Bethe equations. 



4. Bethe Ansatz and Yangian Symmetry 

In this section we will generalize the above construction to arbitrary bound states. We will 
do this by considering coproducts of (Yangian) symmetry generators. This formulation 
allows us to solve (|3.5| ) without knowing the explicit form of the involved bound state S- 
matrix. This will lead to the Bethe equations for arbitrary configurations of bound states. 
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4.1 Single excitations 

We will again start by considering a single excitation in the vacuum 

10) = <)...<). 



(4.1) 



As noted above, it suffices to restrict to two bound state representations. The natural 
generalization of a single excitation wave function is: 



r{K) 



^fc(y) = /(y,Pfc)n'5''"(y'^'')' (4-2) 

Kk 



Restricted to two sites, the wave function is of the form 



(4.3) 



where we again have chosen a particular normalization. The remarkable fact is that one 
can write this as: 



with 



WM := (i^o(j'i,P2)AQi + Ki(pi,p2)AQi) |0), 



(4.4) 



2 1 /\ 1 2 1 2 1 2 1 2 / 



/(P2)5(pi) 



iViPl 



+ 



e 2 



iViPi 



(4.5) 



Xry* ry ' /y» ' 



g^/^ (xg — X]'")(2a;||'"a;^x-|^ — x-^ 



Xi X 



1 -^2 



/(P2)g(pi) 

Vhv{P2) 



f{pi 



iViPi) 



For the moment let us keep /, S arbitrary. The invariance of the S-matrix under Yangian 
symmetry means that 

(4.6) 



In other words, we find: 

S\a) = S (Ko(m,P2)AQi + Ki{pi,p2)Aqi^ |0) 

[Ko{pi,P2)A°Pqi + K,{puP2)A"PQi) S|0) 
[Ko{pi,P2)A°PQi + K,ipuP2)A''Pql^) |0), 
since S|0) = |0). However, we also have 

'Ko{p2,pi)A"PQi + Ki(p2,Pi)A°fQi) |0). 



a 



(4.7) 
(4.8) 

(4.9) 



This means that (^^) corresponds to requiring that Kq and -ftTi are symmetric under 
interchanging pi P2- In other words, to find a well-defined coefficient, we have to solve 



Ko{pi,P2) = Ko{p2,pi), Ki{pi,p2) = Ki{p2,pi), 



(4.10) 
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for the functions / and S. 

It is straightforward to prove that ( [4.1CI| ) is equivalent to the equations: 



Kf{p,) + Gf{p2)S{pi) 
Lf{pi)+Hf{p2)S{pi) 



f{P2) 

f{pi)S{p2), 



(4.11) 



with 



L 



e * 2 



G- 
H 



e 2 



/l2ri{P2) x^-x^ 



(4.12) 



These equations are solved by the f,S found before, i.e. we again find ( 3.11| ) as unique 
solution. The discussion from the previous section holds also here. By fixing the nor- 
malization of the first term to be f{pi), one finds that the second term is factorized and 
fixed. 



Moreover, notice that from this construction we can read off elements of the S-matrix. 
Namely, the coefficients that deal with the scattering of W^^\^ and W^Vw^'^l^ ,. 

We compared the found coefficients with the explicit known S-matrices and we indeed find 
perfect agreement. For example, for they coincide with 09,010,031 and 032, cf. [27]. 

In conclusion, symmetry of the coefficients uniquely fixes the form of our wave func- 
tion. We can now write the wave function, restricted to two sites, completely in terms 
of coproducts and as a consequence ( |3.5|) is automatically satisfied. Finally, the explicit 
expressions for Kq,Ki are 



Koipi,P2,y) 
Ki{pi,P2,y) 




y 



a V ^2 (y-^i ){y-x2) 

4:i I X2 y 

gyxf 



-j^ •^2 2 V I ~'~ 2 ^ ~r "^2 



zjX' y *^2 \ 2 



1 -^2 



X -1 X 



1 -^2 




xt (y - ^1 ){y - ^2 ) 



-J^ ■*-'2 ]^ 2 



2/y» /ya ' nt* ' rv* rv* /y» ' /y^ ' 

1 2 1 2 ■A' ]^ "^2 1 2 



. (4.13) 



This consideration is valid for any bound state numbers and hence wave function (^^) is 
valid for any bound state representations. In particular, all bound state representations 
share the same 5^'^^, and hence that part of the Bethe equations remains the same. 



4.2 Multiple excitations 

When dealing with two excitations, one needs to introduce a level II S-matrix that deals 
with interchanging yi and 7/2- 



Fundamental representations 

Let us first restrict to fundamental representations and reformulate this in terms of co- 
products. The wave function was of the form 

\al3) = \al3)y,y, + ^^^\aP)y,y,, (4.14) 
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where 



S"l«/3)!/i2/2 = M{yi,y2)\aP)y,y,+N{yi,y2)\Pa)y,y,. (4.15) 
The natural way to write this would be: 

{(A,,Qi )(A,,Qi ) + e^pA'y^^y^hl] |0), (4.16) 

with 

A'y^,y^hl := Loiyi,y2,puP2)Ahl + Li{yi,y2,pi,P2)All (4.17) 

By taking a = /?, one easily checks that the first part is indeed of the form (Ay^(Qi)(Ay2Q^). 
Hence, we have to solve Lq, Li such that our ansatz gives 

{f{yi,Pi)f{y2,P2)Siy2,Pi) + Mf{y2,pi)f{yi,P2)Siyi,pi)} W^'^W^^^ 
+Nfiy2,Pi)fiyi,P2)S{yi,pi)W^'^wi^^ 

+e"^/i(yi,y2,Pi)/(y2,Pi)/(yi,Pi)W^i'Vf) (4.18) 

+e'"^h{yi,y2,P2)f{y2,P2)f{yi,P2)S{y2,Pi)S{yi,pi)wi'^W^''\ 

where we keep the functions M, N, h arbitrary. This gives two equations for Lq and two 
equations for Li. Hence, requiring symmetry under pi ^ p2 will give us four equations 
which can be shown to be equivalent to the following set of equations: 

{/l2/2l522 + Mf22hlSl2} = {/ll/22521 + M f 21/12811}^^-^ + N f21 f^Su^^-J^ 

(J 

+ (— /ll/21^121 + /l2/22S'llS'21^122) 

Nf22fuSl2 = {/11/22521 + Mf2lfl2Su}^-^ + Nf2lfl2Su^^ 

(J 

— (— /ll/21^121 + /l2/22S'llS'21^122) ^• 
/ll/2lSl2522/il21 = {/11/22521 + (M - iV)/2l/l25ll}| (4.19) 
+/ll/21^121 ^ 1- /l2/22'S'll'S'2l/il22 ^ 

/12/22/1122 = -{/11/22.521 + (M - iV)/2i/i25ii}| (4.20) 

+/ll/21^121 ^ h /l2/22'S'll-S'2l/il22 ^ , 

where, for convenience, we introduced the short-hand notation fj^i := fiyk,Pi),Ski '■= 
S^^'^{yk,Pl),M := M{yi,y2),N := N{yi,y2) and hijk ■= h{yi,yj,pk). The coefficients 
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B, C, D, E, F are given by 



B 



D 



E 



C = 2irj{pi)r]{p2 



e 2 



X+ (1 - X2 ){XJ 



Xq e 2 



^2 



"16 2 



(4.21) 



e 2 (x-^ (X2 (X;^ — 2x~l) + l)xj' + (xj" + X;^ (x2 x^ — 2))x2 

(1 — X2 — )-'^t 
e~~ (xf -xj")(x^ -X2)(x^ -x+) 

V{P1)V{P2) (1 - x;fX2 - X2 ) 



It is readily seen that these expressions coincide with elements from the fundamental S- 
matrix. Remarkably, these are exactly the same equations that arose in the nested Bethe 
Ansatz. In other words, the coefficients B, C, D, E, F correspond to elements from the 
fundamental S-matrix and we again find ( |3l7| ) as the unique solution for M, A^, h. 

It is worthwhile to note that in this way, we have derived the complete fundamental 
S-matrix. This derivation differs fundamentally from the standard derivations [8, 11] since 
it depends crucially on the full Yangian symmetry. There might be a relation with [44] 
where the fundamental quantum S-matrix was also derived from Yangian symmetry. 

To conclude, let us give the explicit solutions for Lq,Li, 



Ln 



g{yi - y2)x^ x^ 



2i(yi - Xi )(y2 - x^ )(yi - X2 ){y2 - ) 



(4.22) 



rv* T* / T* 1 rf* I r\* I I rf* I 

1 2 1 2 V I 2 1 2 



yiy2Xi X2 



ZjUu^ 1 "^2 



1 2 1 2 ^ *^2 "^1 "^2 1 2 ^ *^2 



l^X]^ -|- X2 X]^ X2 JV'^'l ■''2 -^1X2 ) 1 - "T + "T +1 - I V"''! "''2 "T "''I -^2 



yiy2a;i X2 



(yi - ){y2 - x^ )(yi - X2 )(y2 - X2 ) 



Aig ^x^ X2 x'^X2 



{yi - 2/2) + 



^ju-^ "^2 1 2 



Xi X 



1 -^2 



Note that they are indeed manifestly symmetric under pi <^ P2- 



Bound states 

One might hope that it is possible to repeat the construction of previous section and find 
explicit elements of the bound state S-matrices again. However, when considering bound 

(i) 

states, one encounters a difficulty. There is a new term, which is of the form ^^^11^-2 ■ This 

(i) 

term behaves exactly like W^V^i-i therefore it is hard to rewrite rewrite everything as 
in (4.19) in a unique way and read of S-matrix elements. 
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Nevertheless, we redo the procedure and try to match ( [4.16 ) to the obvious general- 
ization of the two excitation ansatz 



|a/3) = ^MMy2)w[ll . . . w^'X-. ■ ■ ■ • • • + 



(fc) 



.(0 



k<l 



k<l 



+ ^k(.yi)'^kiy2)9{yi,y2,Pk)wX ■■■W^ 

k 



(k) 



Or, restricted to two sites 



(K) 



(4.23) 



(4.24) 



.(2) 



+9{yi,y2,Pi)f{yuPi)f{y2,Pi)W^^^^,^_,wX 
+giyi,y2,P2)fiy2,P2)fiyi,P2)S{yi,pi)S{y2,pi)wXwjX^,^_^ 



.(2) 



(4.25) 



.(2) 



+e"'3%i,y2,Pi)/(y2,Pi)/(yi,Pi)W^2i]i-i^i 
+e'"'h{yi,y2,P2)f{y2,P2)f{yi,P2)Siy2,Pi)S{yi,pi)wXw.^ 



21^2- 



This is indeed possible and imposing symmetry of Lq, Li as before provides equations that 
are uniquely solved by g{yi,y2,Pk) = + M — N) and ( 3.17] ). These factor are, more 

or less, expected from fusion. Plugging these solutions back in Lq,Li, we find the same 
functions Lq, Li as in ( 4.22| ) but one has to bear in mind that now parameterize bound 
state solutions. 

For completeness, we give the explicit two particle wave function (restricted to two 
sites) : 

= {f{yi,Pi)f{y2,P2)S{y2,Pi)+Mf{y2,pi)f{yi,P2)S{yi,pi)}wl\i_,W^%_, 



-Nfiy2,Pi)fiyi,P2)Siyi,pi)W^X.,w;^X-^ 



1 



(1 + M - N)f{yi,p,)f{y2,pi)W^J^,^_,W^ 



2£i 



.(2) 



+^^(1 + M - N)f{y,,p2)fiy2,P2)S{yi,pi)Siy2,Pi)wXw^J^,,-2 



(4.26) 



+e''^Kyi,y2,Pi)f{y2,Pi)f{yi,Pi)W^ll-^W^, 



.(2) 



+e''''h{yi,y2,P2)fiy2,P2)fiyi,P2)Siy2,Pi)S{yi,pi)Wi;;W^;', 

with f,S,h,M,N as in ( ^.11 , 0.17 ). By construction, this wave function satisfies ([j.SD for 
any bound state S-matrix. Hence this solves our two excitation case. In particular one 
finds that also the level II S-matrix, S^^ is unchanged for bound states. 



(l)w(2) 
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4.3 Bethe equations 



By making use of coproducts and Yangian symmetry, we have found a way, independent 
of the explicit form of the S-matrix, to write down Bethe wave functions. This ahowed us 
to find S*^^'^ and we found that the level two S-matrix, S^^ remains unchanged. Hence all 
the higher level factors also remain unchanged. In other words, this yields that the Bethe 
equations for any combination of bound states are given by: 



JPkL _ 



SoiPk,Pl)-'z: 





2 „ X" , 

2 (a) 

nn 

a=l 1=1 




(a) 

/=i Vk 



•^k 



w 



(a) 



+ 



iVk +-(^-^1 



K^i (a) (a) 



n 



W 



(a) 



(a) 

Vk 



1 -L i K™ (a) 



w 



9 

(") _ 2i 

I a 



1 

(a) 
Vk 



g l^k 



W 



(a) 



(4.27) 



with 



^1 _ 1 2i£k X 
H ^k 9 X 



k _ gipfc_ 



(4.28) 



However, note that apart from the parameters x^, the phase factor So{pk,pi) also implicitly 
depends on the considered bound states via [27] : 



-50(^1,^2) 



12 



a{xi,X2) X 



ii-i 



: ^/G{£2-ii)G{i2 + ii) n Gii2 -ii + 21). 



(4.29) 



1=1 



The derived Bethe equations coincide with the equations one expects from a fusion proce- 
dure. This justifies a fusion procedure at the level of the Bethe ansatz equations. 

Finally, we will present the Bethe equations in the mirror theory since they will have 
applications there for the TBA program. In the mirror theory the parameters a, b, c, d 
describing the symmetry algebra have the same dependence on x^ as in the original theory. 
However x^ are now dependent on the mirror momentum 



xf{p) 




-— - i 



(4.30) 



To analyze bound states in the mirror theory it is more convenient to pick up a vacuum 
build up out of fermions (5l(2) sector). The derivation can be repeated for this case and 



-17- 



the Bethe equations are given by [23] 



JPkL 



2 '^(a-) + 





Xll (a) 



n- 

7-1 



(a) 



^+11 (a) , 1 



It; 



(a) 



w 



(a) 



2i 

a 



(a) 



1 

(a) 



If;," + ^ 



(4.31) 



where p is the mirror momentum and are now given by ( [4.30| ). The —1 on the left 
hand side is due to anti-periodic boundary conditions for the fermions [23]. For a detailed 
account of the mirror theory and its analytic properties we refer to [23]. 
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